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\  Abstract 

- 

High  energy  laser  systems  with  highly  accurate  measurements  as  target 
tracking  sensors  use  a  conical  scan  process  to  obtain  a  target  capture  and 
tracking  within  the  narrow  beamwidth.  This  searching  process  and  the  target 
tracking  algorithm  are  major  factors  in  the  performance  of  the  laser  radar/ 
target  tracking  system.  Previous  research  results  use  information-theoretic 
concepts  in  establishing  a  laser  radar/ target  tracking  performance  bound  inde¬ 
pendent  of  the  filtering  algorithm.  A  computer  program  was  developed  to  cal¬ 
culate  the  lower  bound  of  the  estimation  error  due  to  the  non-linear  gaussian 
glint  measurement  process. 

Applying  the  Extended  Kalman  Filter  to  the  angle  estimation  problem  for  a 
gaussian  glint  measurement  process,  the  resulting  filter  is  found  to  have  a 
structure  consisting  of  a  series  of  demodulations  with  gains  adaptively  de¬ 
termined  by  the  resulting  angle  estimate.  The  optimal  performance  of  this 
estimation  process  is  shown  to  be  dependent  on  the  angle  and  a  bound  on  the 
performance  as  well  as  a  numerical  algorithm  is  presented. 


V" 


I.  Introduction 
A.  Mathematical  Model 

The  beam  control  of  a  highly  accurate,  high  gain  laser  radar  system  involves 
an  amplitude  modulafion/glinting/damodulation  in  the  angle  measurement  process. 
The  laser  beam  is  controlled  by  a  mirror  to  minimize  the  angular  error  between 
the  steered  (or  xaser  pointing)  angle  and  the  target  angle  referenced  to  the 
laser  system.  The  laser  radar/target  tracking  system  [1-3]  is  illustrated  .  y 
the  block  diagram  of  Figure  1  (one  dimensional).  The  angular  measurements 
are  developed  by  a  glinting  process  created  by  the  target  and  a  dithered  laser 
beam  described  by  the  following  mathematical  expressions. 

1 .  Dither  Signal 

To  generate  target  tracking  measurements,  the  beam  center  is  dithex-ed  and 
the  returned  signal  is  coherently  correlated  with  that  transmitted.  In  parti¬ 
cular,  the  target/later  beam  angular  difference  with  applied  dithei  signal 
is  modelled  by 


>(t)  =  9b  +  A  sin(utft) 


where 


9b  is  the  target/beam  center  offset  bias,  and 
A  sIn(oj(jt)  is  the  applied  dither  signal. 

2.  Gaussian  Glint 

The  glint  measurement  process  can  be  modeled  by  the  deterministic  gaussian 
function  with  additive  noise 

z(t)  =  g(y,t)  +  n(t) 


where 


g(t)  =  IQ  exp(-y2/2<jg) 


and 


y(t)  is  the  target/beam  center  angular  error, 

Og  is  a  parameter  of  the  glint  return, 

I0  is  the  reflected  intensity  for  zero  error,  and 

-.( t )  is  a  zero  mean,  white  noise  process  with  variance  o“. 

It  is  to  be  noted  that  the  term  "gaussian"  above  refers  to  the  input/output 
process  g(y,t)  and  is  not  stochastic.  The  additive  noise  term,  n(t),  itself 
can  be  a  gaussian  process  but  it  is  not  to  be  confused  with  the  term  "gaussian 
glint"  which  is  the  input/output  functional  relationship. 
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Laser  Radar/Target  Tracking  Stochastic  Control  System 


3.  Signal  Demodulation 


The  output  of  the  gaussian  glint  process  contains  the  information  as  to 
the  magnitude  and  direction  of  the  target/beam  bias.  One  component  of  the 
output  measurements,  z(t),  contains  this  angle  bias,  0^,  as  the  coefficient 
of  the  base  cithered  signal.  Hence,  to  recover  the  bias,  a  simple  demodulation 
process,  commonly  used  in  Communication  Systems  [4]  ,  illustrated  by  Figure  J, 
includes  multiplying  the  output  by  the  dither  signal  and  passing  the  product 
through  a  low  pass  filter. 


z(t) 


sin(uyt) 


Figure  2.  Glint  Bias  Demodulation  Process 


®l 


B.  Statement  of  the  Problem 


The  laser  radar  target  tracking  system  is  a  stochastic  control  system 
whose  objective  is  to  track  a  target  using  noisy  measurements.  The  performance 
criterion  or  measure  of  goodness  of  the  system  is  given  by  both  the  closeness 
of  the  laser  steered  beam  to  the  target  being  tracked  (deterministic  measure) 
and  the  variation  or  jitter  of  the  stei  red  beam  (.stochastic  measure).  In 
many  stochastic  control  systems,  the  optimal  design,  choice  of  parameters,  or 
operating  point  is  usually  a  compromise  between  a  set  of  conditions  which 
yields  the  best  deterministic  response  with  that  which  yields  the  best  stochastic 
estimates.  Although  the  central  problem  addressed  by  this  investigation  is 
that  of  the  stochastic  design,  the  deterministic  nature  of  the  control  system 
must  be  understood  and  considered.  To  this  extent,  the  original  [3]  and  current 
investigation  centers  on  the  problem: 

Given  a  laser/target  tracking  system  with  a  dithered 
signal/gaussian  glint  measurement  process,  evaluate 
the  Information- Theoretic  estimation  performance  bound 
of  the  system. 

The  primary  use  of  an  estimation  performance  bound  is  in  comparison  with: 


1.  the  performance  of  what  is  being  observed  with  an  existing  filter/ 
control  system. 

2.  the  performance  of  an  implementable,  maan-squared-error  (MSS)  optimal 
stochastic  control  which  includes  an  Extended  Kalman  Filter. 
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C.  Proposed  and  Redirected  Investigation 

Reference  [5]  contains  the  initial  proposed  research  investigation  of 
the  estimation  problem  of  the  laser  radar  tracking  system.  Section  II  of  this 
report  contains  the  investigation  results  of  tasks  1,  3,  and  5.  During  a 
meeting  [6]  at  the  Air  Force  Weapons  Laboratory,  preliminary  results  were 
discussed  and  the  investigation  was  redirected  into  the  noise  transfers  both 
in  the  Glinting  and  Demodulation  processes.  Section  III  of  this  report  contains 
the  results  of  this  redirected  investigation. 


There  are  two  key  questions  on  which  this  research  investigation  is  centered:  | 


1.  What  is  the  ultimate  performance  of  the  laser  tracking  system  ? 

2.  What  is  the  performance  of  the  present  process  (is  it  sufficient, 
can  it  be  made  better)  ? 

These  two  areas  were  further  investigated  by  considering  the  extended  Kalman 
Filter  described  in  the  initial  research  investigation  (3]  and  decomposing  it 
into  a  Fourier  Series  representation.  These  results  (Section  IV)  illustrate 
several  fundamental  concepts  regarding  the  present  filtering  process  and  the 
performance  achievable  compared  to  the  ultimate  obtained  through  the  Information- 
Theoretic  approach. 


II.  Information-Theoretic  Approach 


A.  Bounded  Estimation  Performance 


The  initial  investigation  [5]  developed  an  Information- Theoretic  performance 
bound  and  corresponding  numerical  algorithm  used  to  evaluate  an  estimation 
performance  bound  for  the  gaussian  glint  measurement  process  which  will  now  be 
summarized. 

The  estimation  error  entropy  H(x)  of  any  process  is  bounded  by  the  system 
equivocation  H(x|z)  with  conditional  probability  density  function  p(x|z),  i.e., 

~  .  A 

H(x)  >  H(x|z)  =  E(-ln  (p(x|z))) 

which  is  only  a  function  of  the  system  and  independent  of  the  implemented  esti¬ 
mation  process  [7-11]. 

An  important  design  concept  is  that  the  mean-square-error  (MSE)  estimate 
is  a  rainimax  error  entropy  estimate  [10,11],  i.e., 

<v 

min  max  H(x)  <»  min  |  |v(x)|  | 
where  V(x)  is  the  error  covariance  matrix. 

For  the  normal  estimation  problem  [11],  the  derived  optimal  error  entropy  esti¬ 
mation  process  is  the  Kalman  Filter  and  the  error  covariance  matrix  is  identical 
to  the  system  conditional  covariance  matrix,  i.e., 

V(x)  *  V(x | z) . 


I 


i 


-J 
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This  is  very  important  for  understanding  the  Information-Theoretic  approach 
based  on  prior  understanding  of  minimal  variance/Kalman  Filtering  concepts. 

For  the  gaussian  glint  measurement  process  with  the  dither  signal  on  the 
steered  laser  beam,  the  lower  bound  [3,  12,  13]  to  the  error  in  angle  estimation  is 

o2  >  — —  exp(2H( 0))exp(2N1(t))/exp(2N,(t)) 

3  2  ire  1  4 


where  H( 9)  is  the  initial  angle  uncertainty, 


^2  £ 

1  ( t )  =  2a  /  <=xp(-(0b+A  sin(u)T))2/(2aj^+ag))dTi 

an  0,0  0 

2  2  2  2  2V 
°b  =  a0  °g/(2a6  +  ag' * 


N2(t)  =  l 


r=0,2...  n=r  ( r/2) ! ( n-r) ! 2^  0 
and  the  above  integrand  can  be  expressed 


t  1  /  E(g2(i)/a2)kdT, 


t  .  t 

/  {•}  =  [-^  —1  /  exp(-(0b+A  sin(a)T))2/(2o2+o2))dT 


0  a„  °0  n 


where  oc  *  Og/  /k. 

For  the  initial  example,  a  stationary  target  with  parameters: 

A  =  5  normalized  units 
id:  f  =  300  Hz 

Og  =  10 

o9  =  10 

I0/an  s  variable 
0b  :  variable 


For  a  signal-to-noise  ratio  of  unity.  Figures  3a  to  3d  are  the  resulting 
error  variance  lower  bound  for  0b  =  {0,  2.5,  5.0,  10}.  Doubling  the  signal- 
to-noise  results  in  a  faster  responding  lower  bound  as  illustrated  by  Figure  4a 
to  4d,  0b  —  (0,  2.5,  5.0,  10}.  Notice  that  as  the  bias  0b  becomes  much  larger 
than  the  glint  Og,  the  performance  deteriorates. 


B.  Parametric  Studies 

The  following  lower  bound  performance  response  evaluations  were  made  to 
study  the  effects  of  varying  the  system  parameters.  Figure  5  illust *ates  the 
performance  response  as  the  signal-to-noise  ratio  varies  S/N  *  {1.0,  100  and 
10,000}  which  considered  a  unity  dither  amplitude  and  null  bias  angle.  It  is 
obvious  that  increasing  the  signal-to-noise  ratio  improves  the  performance  bound 
response. 


Yorlonc#  Error  Variance 


Figure  5.  Error  Entropy  Estimation  Performance  Bound  for  Signal -to-Noise  Variations 


Figure  6  illustrates  the  case  when  the  dither  amplitude  varies  A  =  {1.0, 
5.0,  and  10.0}  with  a  null  bias  angle  and  a  signal-to-noise  ratio  of  100.  The 
results  illustrated  by  the  figure  shows  that  by  increasing  the  amplitude,  the 
performance  response  is  slightly  slower,  i.e.  approximately  .005  seconds  slower 
for  an  increase  of  10  times  in  dither  amplitude. 

Figure  7  illustrates  the  performance  found  response  as  the  bias  angle 
change  9^  =*  {0.0,  5.0,  10.0}  for  signal-to-noise  ratio  of  100  and  the  dither 
amplitude  unity.  This  illustration  shows  that  as  the  angle  bias  varies  from 
zero  to  10,  the  performance  bound  response  is  approximately  10%  slower. 

In  summary,  the  above  parametric  studies  show  that  the  "best"  (fastest 
performance  bound)  response  occurs  as: 

i)  signal-to-noise  increase 

ii)  bi  a  angle  is  zero 

iii)  small  dither  amplitude 


C.  Information  Flow  Due  To  Demodulation 

For  the  Laser  Tracking  Control  system,  the  demodulation  process  illustrated  by 
Figure  2  consists  of : 

1.  multiplying  the  measurements  z(t)  by  the  fundamental  frequency 
sin(ayt),  and 

2.  passing  the  results  through  a  low  pass  filter  to  obtain  an  estimate 
of  the  tracking  bias. 

It  is  desired  to  determine  if  information  is  lost  during  this  demodulation 
process.  The  measure  of  information  between  the  unknown  bias  0  and  the  measure¬ 
ment  z  is  given  by 

I(z; 0)  -  H(z)  -  H( z | 8) 

where 

z(t)  =  IQ  exp(-y2/2cjg)  +  n(t) 
y(t)  =  9  +  A  sin(n)dt). 

Consider  a  new  measurement 

z | ( t)  »  z(t)  sin(uyt), 

which  Is  the  first  step  in  the  demodulation  estimation  process. 

The  measure  of  information  between  zj  and  9  is  given  by 

I(zi;9)  *  H(zi)  -  H( zj | 9) . 

Using  eq.  (T2.4)  of  reference  {7  or  11],  the  entropy  of  zjCt)  is 
H(z!)  ^  H(z)  +  E{ln| 3z L / 3z | } 

H(zj )  =  H(z)  +  E{ln| sin(uyt) |  }. 


For  the  conditional  uncertainly  term, 

R( zj  |  0)  =»  H(z9)  +  E{ln|sin(uyt)|  |e}. 

Since  the  process  sinCoyt)  is  independent  of  the  unknown  9,  then 
H(zt  |  9)  =  H(zje)  +  E{ln|sin(oyt)  |  }. 

Combining  the  above  two  main  results,  it  is  obvious  that 
I(z.;9)  =  I( z ; 9) 

and  no  information  is  lost  by  the  first  step  of  the  demodulation  process. 

However,  the  second  step  which  involves  low  pass  filtering,  removes  the 
information  contained  in  the  high  frequency  components.  This  point  will  be 
further  explained  when  the  extended  Kalman  Filter  is  decomposed  into  its  Fourier 
Series  components  in  Section  IV.  By  eliminating  the  higher  frequencies  in 
the  measurement  process,  information  about  the  target  is  lost. 


III.  Glint/Filter/Control  System 
A.  Deterministic  Stability 

The  process  of  generating  target  angle  measurements  and  estimates  involves 
a  dithe’’  signal/glint/demodulation  process  is  similar  to  the  amplitude  modula¬ 
tion  communication  technique  [4]  as  shown  earlier  by  Figure  I.  Depending  on 
the  operating  or  glint  bias  point,  the  received  signal  contains  the  information 
on  the  direction  and  magnitude  of  the  angle  pointing  bias.  The  nature  of  the 
filter/control/mirror  deterministic  portion  of  the  feedback  system  can  be 
analyzed  using  classical  techniques.  The  modulation/glint/demodulation  process 
can  be  linearized  using  the  first  order  approximations  of  the  glint  output  with 
respect  to  the  angle  bias  (neglecting  the  dither  signal  carrier)  as  shown  in 
Figure  8 


where 


Figure  8.  Linearization  Model  of  the  Gaussian  Glint  Process 


Effective  Gain:  K0  *  -  G0  9j,/Og 


DC  Level: 


G0  3  T-o  exp(-'iJ/20^) 


Average  Bias:  9*, 

(operating  point) 
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Another  method  to  determine  the  '•ffective  system  gain  of  the  glinting  process 
is  to  use  the  Fundamental  Describing  Function  [14].  In  either  cat,e,  it  is 
obviously  a  function  of  the  bias  operating  point.  The  standard  modulation 
technique  consists  of  using  a  sinusoidal  signal  to  carry  the  information,  and 
in  our  case,  it  contains  bias  error.  The  bias  error  can  be  obtained  by  a 
coherent  detection  process  which  consists  of  multiplying  the  glint  output 
(measurement)  bv  the  same  sinusoidal  signal  and  sending  the  resultant  signal 
through  a  low  pass  filter.  This  process  is  essentially  the  mathemr.uical  process 
used  to  determine  the  basic  Fourier  Series  harmonic  in  a  periodic  signal. 

This  concept  will  be  used  in  Section  IV  when  the  Extended  Kalman  Filter  is  decom¬ 
posed  into  its  basic  components.  The  overall  effective  system  gain  will  be 
maximum  when  the  signal  is  at  the  gaussian  glint  point  <jg  and  vanished  at 
zero,  hence,  the  system  gain  varies  according  to  the  bias  point. 

Table  1  lists  the  various  transfer  functions  [6,15]  of  the  feedback  path. 

The  transfer  function  H(jui)  of  Figure  9  illustrates  tne  general  nature  of 
the  filter/control/mirror  path 

Gf (s)  =  F(s)  C(s)  M(s) 

F(s)  =  G2(s)  G^s)  (without  the  notch  filter  G3(s)) 

C(s)  *  1/s  (without  compensator) 

The  effect  of  the  notch  filter  G^(s)  adds  zeros  at  the  dither  and  twice  the 
dither  frequency  (300  and  600  Hz)  and  is  intended  to  remove  frequency  modes 
which  are  caused  by  non-linearities  and/or  a  non-perfect  demodulation  process. 
However,  by  observing  Figure  9,  the  feedback  system  G(s)  is  basically  a  Type 
II  low  pass  filter  with  an  essential  triple  break  at  150  Hz.  The  suppression 
effect  of  this  filter  on  a  signal  at  the  dither  frequency  is  large  enough  so 
that,  for  analysis  purposes,  the  notch  filter  Gyis)  will  not  be  included.  For 
example,  there  is  an  approximate  30  dB  transfer  function  decreased  froma  fre¬ 
quency  increase  of  150  tc  300  Hz  which  represents  a  1/25  relative  signal  change. 

The  effect  of  the  compensator  Gc(s)  essentially  stabilizes  the  system. 
Figures  10a  and  10b  show  that  Bode  Plots  (amplitude  and  phase)  of  the  feedback 
system  with  and  without  the  compensator  which  is  used  to  provide  positive  gain 
and  phase  margins  (gain  margin  of  75  dB  at  35  Hz). 


3.  Stochastic  Analysis 

The  demodulation  process  consists  of  multiplying  the  measurements  with  the 
initial  dithered  signal  and  passing  the  product  through  a  low-pass  system  which 
consist?  of  a  Filter/Control  subsystem.  In  particular,  consider  the  process  of 
the  demodulation  process  as  shown  in  Figure  2  where  the  noisy  measurements  is 
multiplied  by  sin(<iyt)  and  then  passed  through  a  low  pass  filter  to  regain  the 
D.C.  level  which  contains  the  glint  bias.  The  noise  analysis  of  this  demodula¬ 
tion  process  co  .siders  the  input  noise  n^(t)  applied  to  the  demodulation 
process  is  characterized  by  its  spectral  density  of  Figure  11. 
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Figure  11.  Spectral  Density  of  the  Input  Noise  to  the  Demodulation  Process 


The  autocorrelation  of  this  spectrum  is 


Rj^Cx) 


(2irr)‘ 


1  -  cos(u)ix)  COS(u)ot)  -  COS((D9T) 

{  - i  * - 1 - }  . 


f3  -  f2 


Its  variance  is  R^(0)  «  A[f3  +  f2  ~  fj]  which  is  the  same  value  obtained  when 
the  spectrum  Si(f)  is  integrated  over  all  f 

00 

*  \(0)  =■  /  SjU)  df. 

i  —co 

Now,  the  output  noise  process  can  be  modeled  by 
n0(t)  =  n^t)  sin(ayt  +  6) 

where  n*  has  Ri(x),  Si(f)  general  characteristics,  in  particular,  that 
specified  above.  Since  the  input  noise  is  purely  random  in  nature,  it  is 
uncorrelated  with  the  modulation  term  sin(uyt)  or  when  it  begins.  Therefore, 
we  can  assume  that  the  modulation  term  has  an  arbitrary  phase  which  is  uniformly 
distributed  between  0  and  2 it.  The  output  correlation  function  is 

Ro( t)  =  E[n0( t)  n0( t+x) ] 

=  E[n^(t)  sinCuyt+'Ji)  ni(t+x)  sin(uy  (t+x)+ij>)  ] 

*  E[n^(t)  n^(t+x)]  E[sin(oyt+i|>)  3in(uy(t+x)+<p)] 

■  Ri(x)  E[cos(uj<4x)  -  cos(uid(2t+x)+2i!i)]/2 

Rq( x)  *  Ri<x)  zos(u4t)/2. 

The  spectral  density  of  the  output  Sc(f)  is  the  Fourier  Transform  of  R0(x>. 


s0(f)  =  1/2  /  R±(x)  co,'>Udx)  dx 


...  -j(or-wd)7  -j  ( or*-u>d)  T 

=  1/4  /  R^( x)  [e  +  e  ]dx 

—  00 

s0(f)  =  [Si(f-fd)  +  Si(f+fd)]/4. 

Hence,  the  output  spectral  density,  S0(f),  is  1/4  the  sum  of  two  densities; 
one  of  which  is  the  input  density  shifted  left  fd,  and  one  shifted  to  the 
right  £d 

Depending  on  the  value  of  fd  relative  to  f^,  f£  and  f 3 ;  the  shape  of  the  output 
spectrum  will  have  different  characteristics.  For  example,  Figure  12a  shows  th< 
output  spectrum,  S0(f)  for  f  +f 

-V-  '  fi  <  f2 


Figure  12a.  Spectral  Density  of  the  Noise  After  Multiplication 
With  the  Dithered  Signal  -  Large  f<j 


f A  *  *2-*d 


fB  -  f3-fd 

*c  -  Vfi 


fF  ’  fl+fd 
f F  ■  f2+fd 


s0(f)  =*  1/2  /  Rd(x)  ccK'.(u)dT)  e"^1  dx 


/  r*  ,  v  r  ~j ( tAj-aid) T  -j(uH>d)T 
1/4  /  R^x)  [e  +  e  ]dx 


SQ(f)  =  [S±(f-fd)  +  Si(f+fd)]/4. 

Hence,  Che  output  spectral  density,  S0(f),  is  1/4  the  sum  of  two  densities; 
one  of  which  is  the  input  density  shifted  left  fd,  and  one  shifted  to  the 
right  fd 

Depending  on  the  value  of  fd  relative  to  f^,  f£  a°d  f 3 ;  the  shape  of  the  output 
spectrum  will  have  different  characteristics.  For  example,  Figure  12a  shows  the 
output  spectrum,  S0(f)  for  fi+f-> 

J_i  <  fd  <  f2 


SoCO 


Figure  12a.  Spectral  Density  of  the  Noise  After  Multiplication 
With  the  Dithered  Signal  -  Large  fd 


f3+fd1 


fA  "  f2*fd 
f B  *  f 3-f d 
fC  "  fd“fl 


f R  ’  fl+fd 

fp  =*  f2+f d 
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The  variance  of  the  output  is  the  area  under  S0(f)  which  is 

°out  “  2{  N[ f 2  -  fd3/2  -  Kfj/4  +  Nff3  +  fd  -  (f2  -  fd)]/4} 

«  N[f 2  -  fd  -  f i/2  +  (f3  -  $2  +  2fd)/2] 

-  Nl  f 3  +  f2  -  f  j  ]  /  2 
°out  =  °in^2« 

The  output  variance  is  1/2  the  input  variance,  but  its  spectrum  has  shifted 
with  a  flat  spectrum  about  f=0  whereas  the  input  has  a  "V"  about  f=0. 

fij.f. 

Figure  12b.  shows  the  output  spectrum,  SQ(f),  for  <  - - 


Figure  12b.  Spectral  Density  of  the  Output  Noise  After  Multiplication  With 
the  Dithered  Signal  -  Small  fd 

fA  "  fd~fl  fE  =  f3*fd 

fB  *  fl+fd  fF  =  f2+fd 

*C  ' 

The  variance  of  the  output  is-,  again,  the  area  under  S0(£): 

■&c  -  2  f  (<f2+fd>  +  (f2-fr2£d>  +  <fd-f.5  +  <f3-f2>  +  Ell 


_[f3+f2-f j] 
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i 

j 


cout  “  0in^ 


Again,  Che  spectrum  is  flat  about  f=0  whereas  the  input  has  a  "V"  about  f=0. 

This  output  spectrum  is  then  passed  through  a  low  pass  filter.  The  noise 
spectrum  of  the  feedback  output  is 

SfQ(f)  =  |H(jm)|2S0(f) 

where  H(jaj)  is  the  combined  transfer  function  of  the  f liter/ control/mirror  feed¬ 
back  subsystem.  There  are  two  types  of  feedback  systems,  Type  I  and  Type  II 
systems  to  be  considered.  With  a  Type  II  system, 


Hn(s) 


KH0(s) 

s^ 


where  H0(s)  is  the  Type  0  portion  of  the  subsystems.  Assuming  that  K 
and  since  the  effects  of  H0(s)  is  negligible  at  low  frequencies,  then 
spectrum  of  the  output  filter/control  is 


stc(f) 


S0<*>  # 

(2ir)4f4  * 


=  1 
the 


Approximating  SQ(f)  to  be  uniform  between  ±  f2,  then 


Sfo<f> 


N 


2(2TT)4f4 


5  *f.  ^  f  ^  f' 


The  autocorrelation  function  of  the  feedback  output  is  the  inverse  Fourier 
transform  which  is 


Rfo(  t)  =  l/2ir  Afo(f)  dm 
—  00 


o>2 

=  1  / 2 ir  /  N/2ui4  du>  . 

—  u>2 


The  noise  variance  at 
the  spectrum 


°fo 


the  feedback  output  is 


jf2  N  df 
0  2(  2ir) 4  f4 


Rfo(0)which  is  the  area  under 


*2 
0 

and  a  limiting  process  must  be  applied 


-N 

3(2r)4f3 


- £—  [lim  1/a3  -  1/f3) 

3(2ir)^  a+o 


the  above  becomes  unbounded  as  a  -*•  0.  This  ioplies  that  the  realization  of 
the  noise  may  be  any  value. 

Consider  the  case  where  we  have  a  type  I  system 
Hj(s)  =  KH0(s)/s. 

As  before,  assume  K  =  1  and  since  the  effects  of  H0(s)  is  negligible  at  low 
frequencies,  then  the  output  spectrum  is 


S.(f) 

Sf0(f>  = 


(  2  rr) 2  f  2 

Again,  as  before,  if  S0(f)  is  uniform  in  the  region  -f2  <  f  <  f2*  then 


sfo(f>  • 


N 


>  — f  j  f  <  f  2  • 


2(2ir)2f2 

The  noise  variance  is  the  area  under  this  spectrum 

.f  2 


N 


df 


°f°  '  2  !°  2<2«)2  £2 


-N 


(2  it)  2  f 


0 


Again,  as  before,  the  above  is  evaluated  using  a  limiting  process; 

N 


4 


(2ir)<-  a-K> 


{lim  1/a  -  1 /f 2 } 


which  becomes  unbounded  and  the  same  results  hold. 


C.  Signal  and  Noise  Transfer  for  the  Gaussian  Glint  Process  (about  y  =  0,  og) • 

The  signal  and  noise  transfer  characteristics  of  the  gaussian  glint  process  can 
be  seen  by  the  following  approximation  process.  Consider  a  Taylor  series  ex¬ 
pansion  of  the  gaussian  glint  process  about  some  point  yQ  =  0  due  to  a  small 
signal  and  noise 


y(t)  =  0  +  A  sin(uyt) 


00  i 

f (y)  =  exp(-y2/2a2)  =*  £  — 

8  i-0  ay1 


(y  -  9)1 

i! 


y«9 
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The  various  (first  4)  derivatives  of  the  function  f(y)  are 


f (y)  =  exp(-y2/2<^) 
f ' ( y )  =  (-y/a2)exp(-y2/2a2) 
f”(y)  =  (y2  -  a2)exp(-y2/2o2)/og 
f'"(y)  =  (3ya2  -  y3)exp(-y2/2a2)/Og 
f”"(y)  =  (y4  -  6y2o2  +  3a4)exp(-y2/2ag)/cjg 
There  are  two  operating  points  of  special  interest 


zero  slope  f'=  0  ==>  y  =  0 
maximum  slope  f"  =  0  ==>  y  =  ±0g 


Evaluation  of  the  various  derivatives  under  zero  and  maximum  slope 
y  0  ;  f(y)  =  e°  =  1 
f’(y)  =  0 


f"(y)  =  -  l/o: 


g 


f"’(y)  =  0 


f” "(y)  =  3/ a4 


f(y)  =  l//e 
f’(y) 
f"(y)  =  0 


-  l/ag  /e 


f"' (y)  =  2/cl  /e 


f""(y) 


-  *>4 


In  Tabular  form 


y  =  0 

y  =  ag 

f(y) 

1 

l//e 

f '  (y) 

0 

-1/og/e 

f"(y) 

i  '1/og 

0 

f” (y)  | 

1  ° 

f" "(y) 

3  *4 

-2/cUe 

Maximum  Slope: 

Consider  the  maximum  slope  bias  point.  Using  the  first  4  terras  of  the  series 
expansion  about  y0  =  ag  with  y  =  ag  +  A  sin(uyt)  +  tHU)  where 
ni(t)  is  the  input  noise,  the  output  is 

z(t)  =  l//e  {l  -  (A  sin(ait)  +  ni(t))/og  +  0/2! 

+  2(A  sin(ajt)  +  ni(t))^/Og  3!  +  ••••  } 


The  most  significant  terras  of  the  expansion  are: 
signal  {l  -  A  sin( wt)/cg }/ /e 


noise  n^(t)/0g/t- 


Hence  both  the  signal  and  noise  process  are  preserved  with  a  change  in  amplitude. 
It  the  maximum  slope,  the  expected  signal  and  noise  terms  are  easy  to  understand 
and  observe  using  Ihe  above  simple  example  as  illustrated  by  the  frequency  spectrum 

of  Figures  13a  and  13b. 


Minimum  Slope: 

Using  the  first  4  terms  of  the  series  about  0  =  0  the  output  is 

Z(t)  =  1  +  0/1!  -  (A  sinUt)  +  n^/o2  2!  +  0/3!  +  3(A  sin(wt)  +  4!  +  ••• 

2  2 

=  1  -  [A2sin2(ujt)  +  2A  sin(uJt)  n^(t)  +  n^  (t)]/2og 
+  [A^sin^(arr)  +  ...  +  n^^(t)]/80g  +  ••• 


The  most  significant  terms  of  the  output  are: 

signal  1  -  A2  sin2(wt)  /2a2 

noise  [2  n^t)  A  sin(  >t)  +  ni2(t)]/2o2 

Hence  both  the  signal  and  noise  process  change  at  the  zero  glint  bias  point. 

IZ  to  the  extremeSnon-linearity  at  this  point,  both  the  output  signal  and  noise 
change  characteristics  the  modulation  signal  amplitude  is  now  at  twice  the 
original  frequency.  T'ne  output  noise  has  two  components,  one  modulated  by  the 
du£  frSS  a^d  one  .odulat ,d  by  itself.  These  wo  modulated  noise  terms 
will  give  rise  to  spectrum  components  not  present  at  the  input  as  illustrated 
by  Figure  13c. 
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Figure  13a.  Frequency  Spectrum  of  a  Signal  and  Noise  at  the  Input  of  the 
Gaussian  Glint 


-*d  fd 

Figure  13b.  Frequency  Spectrum  of  the  Output  of  the  Gaussian  Glint  Process 
with  the  Operating  Point  0^  =  og 


2fd 


Frequency  Spectrum  of  the  Output  of  the  Gaussian  Glint  Process 
the  Operating  Point  6^  *  0 


IV.  Extended  Kalman  Filter  for  Laser  Tracking 
A.  Extended  Kalman  Filter 

A  straight-forward  application  of  an  Extended  Kalman  Filter  [3j  for  the 
bias  angle  estimation  yields  the  following  set  of  equations  to  be  implemented. 

0b  =  K(t) ( z( t)  -  g(y,t)l 

y(t)  =  6b  +  A  sinCujt) 

g(y,t)  =  IQ  exp(-y2/2t^) 

K(t)  =  -  !|X_  g(y,t) 
an  °g 

=.  -  V-|  4  g2(y,t) 

°n  °g 

Note  that: 

h  *  -  K2(t)°n 

Figure  14.  is  a  block  diagram  for  the  above  filtering  process. 


B.  Fourier  Series  Decomposition  of  the  Extended  Kalman  Filter 

The  following  will  be  a  decomposition  of  this  filter  into  the  Fourier 
Series  representation  of  the  filtering  process.  As  shown  earlier,  the  Taylor 
Series  expansion  of  the  gaussian  glint  process  about  a  point  0  is 


«<y>  -  l  & 

i=*0  ay1 


y=  0 


i! 


With  the  input 


y  =  0  +  A  sinCoyt), 


then 


-0Z/2q2 


62  '  oi 


f(y)  =  e  g  {l  -  -2  A  sin(codt)  +  — - — A2sin2(mdt) 


30a2  -  03 


2! 

04  -  6  02a2  +  3a4 


°g  3! 


A3sin2(a)dt)  +  - g-  —£ - S.  A^sin^Cw^t)  -1  ...  } 


°g4! 


Figure  14.  The  Extended  Kalman  Filter  for  the  Gaussian  Glint  Process 


The  Kalman  filtering  gain  is 
-V~  y 

K(t)  — § —  g(y) 

<*n  o| 

where  Vg  is  the  error  variance 

g(y)  =  I0f(y),  and  y  =  e  +  A  sinC^t), 


-V~  -02/2o|  .  , 

K(t)  =  — § -  [0  +  A  sin(a)dt)]  IQ  e  *  l*}  * 

On  Og 


_V~i  e"92/2a2 
L(V~,0)  =  — s 


,  then 


°n  °1 


^2  2 

K(t)  =  L  [0  +  A  sin(u)dt)]  {l  -  — 2  A  sin(o,dt)  +  — A2sin2(o>dt) 


3 0g2  _  @3  _  01*  -  662a~  +  2ag  ,4.^4^  t-\+  l 

+  i!fS _ _  A3sin3((odt)  +  - r- f - *  A  sin  (wdt)+*“f 

4  31  as  4- 

L  {CQ  +  CL  sinUdt)  +  C2  sin2vu>dO  +  C3  sin3(a)dt)  +  C4  sin4Udt)  + 


where 


Cn  -  0 


ol  -  02 


C1  -  A  [  -5-2 -  J 

1  4 

.2  93  -  0a|  2 ! 0 

C„  =  A-  [  — 5 — S-  -  -Y-  ] 

2  2!  °g 

2  03  -  30a2 

-  A-  [  _  ■  8  1 

2!  ^ 


C,  -  ^  l 


*3  302a2  -  04 


02  -  a2 


a3  -8-  +  68*o|  -  3<g 

_  * '  *  K 


O 

,4  05  -  603a2  +  30Cp  4!  30a2  -  03 

r  =  **  r  ,  P  --**  + - 8 — ? - 

c4  rr 1  «§  3!  oi 


.4  05  -  lO03a2  +  I3  0a£ 

A  f  _ _g  - s.  J 
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dsing  the  various  identities  for  sinn(a),  i.e. 

sin1 2 *(a)  =  [1  -  cos(2a)]/2 

sin2(a)  -  [3sin(a)  -  sin(3a)]/4 

sin4(a)  =  [2  -  2cos(2a)  +  icos(  4ct)  ] /4  , 

2  2 

the  first  few  terms  of  tae  Fourier  series  representation  of  the  Kalman  gain  are 
4  4 

K(t)  =  Lfa0  +  l  cos(i<udt)  +  £  sin(iu^t)! 
i=l  i=l 

ao  =  Co  +  C2  +  }  C4  +  * '  * 

aj  =  0 

a~  =  -  —  C0  -  —  C/  +  . . . 

2  2  2  2  4 

a3  =  0 

a4  =  |  c4  +  *  ‘  * 

bl  =  C1  +  1  C3  +  ... 

b2  =  0 

b-j  =  -  —  +  ... 


b4  =  0 

The  Kalman  gain  is  used  to  develop  the  time  derivative  of  the  tias  estimate  by 
multiplying  it  with  the  measurement  residual 


z(t)  =  z(t)  -  z(t) 


0(t)  =  K(t)  z ( t ) 

Figure  15a  illustrates  the  process  of  obtaining  the  estimate  of  the  glint  bias 
and  Figure  15b  is  the  process  now  being  used  to  generate  the  estimate.  It  is 
obvious  that  the  present  process  is  a  simplification  of  the  optimum  process  in 
that: 

1.  only  the  information  in  the  fundamental  harmonic  is  used,  and 

2.  the  transfer  gain  is  constant  instead  of  the  time  varying 

optimum  process. 
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Kalman  Process  in  Fourier  Series  Form  to  Estimate  the  Glint  Bias 


The  second  term  of  ..he  Kalman  Process  illustrated  on  Figure  15a  is  similar  to 
the  present  process  of  Figure  15b.  In  fact,  the  gain  is  an  adaptive  process 
of  the  glint  bias 


L  b. 


~  -0?72o| 

V9  IQe  s 


«f 


<C1  +IC3> 


Normalizing  this  bias  term  using  d  =  0/ <jg 
Vo  Ine-d2/2 

[  A(  1  -  d  )  + 

8o| 


Lb,’ 


'8  Lp* 
°n  Og 


t A(  1  -  d2)  +  (-  d4  +  6d2 


-  3)] 


V~  I 

liA.  f(d) 

On  o| 


where  the  following  section  describes  the  normalized  gain  function  f(d)/og 


C.  Coefficient  Evaluations  for  the  Fourier  Series  {a^,  b^} 

The  Fourier  Series  Coefficients  {a^,  bjJ  for  the  Kalman  gain  can  be 
evaluated  using  a  glint  bias  to  standard  deviation  normalized  variable 


d  *  8/og 


This  process  will  show  the  relative  importance  of  each  coefficient.  The  various 
coefficients  are 

a  =■  da„  +  -  3d)  +  (d5  -  10d3  +  15d) 

°  g  4og  64  o| 

a^  *  0 

a~  =  -  —  (d3  -  3d)  -  (d5  -  10d3  +  15d) 

2  4°g  48  °! 

a.^  ~  0 

a,  =  Jl.  (d5  -  10d3  +  15d) 

4  192 o| 

b ,  =  A(  1  -  d2)  +  d4  +  6d2  -  3) 

^2  -  0 

b,  =  -  d4  +  6d2  -  3) 

3  24°! 

=>  0 


Figure  16a  and  16b  lists  normalized  values  for  these  coefficients  for  various  * 
operating  points  d  ■  {-1  to  1}  which  corresponds  to  a  variecy  of  glint 
bias  8b  points  and  for  values  of  dither  amplitude  A  =  l/2ag,  and  <jg. 
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Normalized  Angular  Bias 


D.  Optimal  Estimation  Error  Variance  (Bound) 

The  following  will  be  a  bounded  closed  form  solution  for  the  error  variance. 
The  error  variance  (using  V  for  V~)  differential  equation  is 

6 


2  2  2 

v  =.  _  V  y  g  (y,t) 

°n  °g 

where 

y(t)  =>  0  +  A  sin(ajdt),  and 

g(y,t)  =  I0  exp(-y2/2a|) , 

then 

•  2  2 

V  =  -cV  [0  +  A  sin(o)dt>] 

exp(-y2/ Og) • 

where 

I2 

c  *  0  (constant). 

On  Og 

Separating 

variables ,  the  above  can  be 

expressed  as 

*  cy 2  ( t )  exp(  -y 2  /  Og )  dt . 


Upon  integrating, 

-7—  -  —  *  c /  y2(T)exp(-y2(T)/Op)dT 
V(t)  VQ  '  g 

where  V0  is  the  initial  condition  on  V.  Since  the  exponential  function  can  be 
bounded 

exp[-( | 0]  +  j A | ) 2 / ]  <  exp(-y2(t)/og)  <  1, 

the  above  integral  can  be  bounded  by 

exp[-lii|i  +.  JAB,  .]  FCc)  <  c  /  y2(T)exp(-y2/o2)dT  <  F(t), 

1 

where  F(t)  is  a  monotonically  increasing  function  of  time 
t  9 

F(t)  -  c  /  [e  +  A  sin(uidT)]zdT 
o 


2  2  2 
c  {(0  +  A  / 2 ) t  +  2A0(l-cos(a),jt))/a)d-A  sin(2u){jt)/4a)d}, 


The  error  variance  can  now  be  bounded 
V 


_ _  <  V(t)  <  _ 2 _ 

1  +  VoF(t)  1  +  V  F(t)exp[-( I 0|+ j A| ) 2/o2] 
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2.  Optimal  Estimation  Error  Variance  (Numerical  Evaluation) 

1.  Numerical  Solution 

The  differential  equation  of  error  variance  for  the  Extended  Kalman  Filter 
[4]  can  be  expressed  as 


where 


and 


V  =  -  c  V2  y2(t)  exp(-y2(t)/o|) 


c 


t2 

*o 


a 


4 

g 


y(t)  =  0fc  +  A  sir(ayt). 

Separating  variables  and  integrating  yields 

-  —  =  c  y2(t)  exp(-y2(t)/a|)dt 
v2  8 

-L-  =  c  /  y2(  r)  exp(-y2(  t)/  o|)  dT  + 

V(t)  o  g  V0 

and  upon  simplification,  the  error  variance  become 

V(t)  =  - - 

1  +  c  VQ  /£{•} 

where  V0  is  the  initial  error  variance.  Since  the  integrand  is 


(•}  •  y2(t)  exp(-y2(t)/ a|) 

is  a  periodic  function,  then  the  integral  can  be  evaluated  at  specific  times 
which  are  integer  multiples  of  ,  the  period  of  the  dithered  signal. 


{•}  =  n 


I' 

o 


{•}  • 
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The  above  one-period  numerical  integration  process  is  easily  obtained  with  Oe 
use  of  a  sufficiently  small  integration  step  size  (At). 


2.  Parametric  Studies 

Similar  to  the  performance  studies  from  the  Information-Theoretic  Bound 
(part  II. B.),  Figure  17  illustrates  the  numerical  evaluation  of  the  error 
variance,  V(t),  for  signal-to-noise  ratios  of  S/N  =  {1.0,  100,  and  10,000} 
for  a  unity  dithered  amplitude  and  a  null  bias  angle.  The  strong,  direct  in¬ 
fluence  of  the  S/N  is  apparent  from  the  response  illustrations  of  Figure  17 
and  the  integral  solution  of  V(t)  where  the  multiplicative  parameter  "c" 
contains  the  signal-to-noise  ratio. 

Figure  13  illustrates  the  error  variance  response  for  various  levels  of 
dither  amplitude  with  S/N  =  100  and  9^  =>  0.  The  response,  V(t),  is  shown  to 
improve  (decreases  faster)  as  the  dither  amplitude  increases,  however,  the 
rate  of  improvement  decreases  when  at  large  values  (i.e.  from  A  =  5  to  A  =  10). 

Varying  the  bias  operating  point,  Figure  19  illustrates  an  improvement  in 
the  response  V(t)  as  0^  increases  with  a  limiting  effect  at  higher  levels  of 
bias  (i.e.  from  0^  *  5  to  9^  =  10). 


3.  Comparison  of  the  Kalman  Error  Variance  and  the  Information-Theoretic  Bound 

Two  approaches  have  been  used  to  evaluate  the  performance  of  the  optimal 
angle  estimation  for  the  gaussian  glint  process.  Although  the  Information- 
Theoretic  and  the  Extended  Kalman  results  are  related  to  the  same  quantity, 
their  Interpertation  and  the  assumptions  made  for  their  evaluations  are  dif¬ 
ferent.  However,  they  both  give  some  measure  as  to  the  overall  performance  of 
the  optimal  angle  estimation  process  as  illustrated  by  Figure  20  which  is  noth¬ 
ing  more  than  combining  Figures  5  and  17.  The  results  illustrate  the  lower 
bound  attribute  of  the  Information-Theoretic  evaluation  and  the  expected  re¬ 
sponse  when  the  Extended  Kalman  Filter  is  used. 


-37- 


Time  (Seconds) 


Figure  18.  Kalman  Error  Variance  for  Dither  Amplitude 


Kalman  Error  Variance  for  Bias 
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Comparison  of  the  Informat  ion- Theoretic 


V. 


Conclusion 


A.  Summary  of  Results 

An  Information-Theoretic  Bound  was  found  for  the  Laser  Radar  Tracking 
Problem.  It  was  shown  that  the  signal-to-noise  ratio  strongly  affects  the 
performance.  There  is  also  a  slight  variation  in  performance  as  a  function 
of  angle  bias  and  signal  amplitude  but  not  a  significant  one. 

An  Informatic  .-Theoretical  Analysis  of  the  demodulation  process  showed 
that  there  was  no  loss  of  information  due  to  the  first  step  of  the  process, 
multiplication  of  the  gaussian  glint  returned  by  the  transmitted  dithered  signal. 
However,  the  suppression  of  the  higher  frequency  components  by  the  second  step 
of  the  demodulation  process,  low  pass  filtering,  destroys  information  contained 
the  glint  return. 

The  type  II  feedback  transfer  function  is  shown  to  be  stabilized  by  the 
compensator  where  a  deterministic  positive  gain  margin  is  given.  The  first 
step  of  the  demodulation  process  was  shown  to  redistribute  the  noise  spectrum 
contained  in  the  output  gaussian  glint  to  significant  levels  at  the  low  fre¬ 
quencies  (D.C.  region).  Coupling  this  with  the  type  II  (or  type  I)  feedback 
transfer  function  causes  an  unbounded  stochastic  variation. 

The  bias  operating  point  has  a  strong  effect  on  the  spectral  transfer 
function  of  the  non-linear  gaussian  glint  process.  In  particular,  a  small  signal 
analysis  showed  that  the  zero  bias  operating  point  changes  the  deterministic 
and  stochastic  spectral  response  significantly.  When  the  operating  point  is  at 
the  maximum  slope,  the  output  spectrum  reflects  the  input  spectrum. 

The  Extended  Kalman  Filter  was  applied  to  the  open  loop  angle  estimation 
problem  involving  a  gaussian  glint  measurement  with  a  dithered  applied  signal. 
Decomposing  the  filter  formulation  into  a  Fourier  Series  representation,  it  was 
shown : 

i)  the  angle  estimation  process  consists  of  a  series  of  demodula¬ 
tion  processes  corresponding  to  various  frequencies  of  the 
applied  dither, 

ii)  the  gains  are  chosen  adaptively  according  to  the  bias  angle,  and 

iii)  the  overall  optimal  error  variance  decreases  monotonically  as 
time  increases. 

The  optimal  adaptive  filtering  process  presented  is  in  the  form  similar  to 
the  present  system  and  improvements  can  be  made  by  using  additional  elements. 
Bounds  on  the  optimal  error  variance  were  established  in  closed  form  and  a 
numerical  technique  was  also  presented. 

The  error  variance  for  the  angle  estimation  is  shown  to  be  bounded  and  a 
numeric  process  is  also  presented.  Evaluating  the  error  variance  for  various 
parametric  cases  show  a  strong  dependence  on  signal-to-noise.  It  is  also  shown 
that  ;he  dither  amplitude  and  bias  operating  point  has  an  effect  on  the  estima¬ 
tion  performance. 
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B.  Extension 

There  are  a  number  of  natural  extensions  to  the  material  presented  in  this 
paper : 

i)  Spectrum  Analysis  of  the  System 


Items  in  this  investigation  would  include  both  signal  and  noise 
transfers  through  the  non-linear  elements  of  the  gaussian  glint 
and  the  demodulation  processes.  The  spectrum  analysis  of  the 
closed  loop  controller  and  the  Extended  Kalman  Filter  will  also 
be  included. 

Parameter  Variations  of  the  Stochastic  Control 

There  are  various  operating  points  and  system  parameters  that 
yield  optimum  performance.  A  sensitivity  study  of  these  design 
parameters  and  operating  points  will  be  included. 

Adaptive  Filter  Implementation 

There  are  two  forms  of  the  Extended  Kalman  Filter  to  be  imple¬ 
mented.  Each  of  these  forms  should  be  investigated  as  to  their 
complexity  of  implementation.  Also  included  in  this  investigation 
should  be  a  sub-optimal  filter  which  consists  of  the  main  element 
of  the  Optimal  Filter.  A  performance  analysis  of  this  design  will 
indicate  the  effects  of  the  sub-optimal  design. 

Optimal  Controlling  Process 


The  control  process  of  the  present  system  consists  of  a  simple 
integrator  on  the  bias  error.  An  optimal  design  based  on  modern 
techniques  which  includes  performance  criterion  of  minimum 
settling  time,  control  energy,  jitter,  etc.  would  be  included. 
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